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Abstract 

Electric potentials for two-dimensional cross-sections of ion focusing 
devices used in a mass spectrometer are calculated via a series of computer 
programs designed to compute potentials between areas of fixed voltages. Ion 
trajectories within these devices may be determined by computer and a histo- 
gram obtained which plots ion density against ion position along a plate of 
the focusing device. For each lens system, a plate voltage may be changed, 
the electric potentials recalculated, and a new histogram calculated in or:ler 
to determine if the new voltage configuration has increased the device's effi- 
ciency. This process may be repeated until the optimum voltage values have 
been found for maximum particle transmission in each focusing device. 

I. Introduction 

A sequence of computer programs has been written by Dan Anderson which are 
capable of determining electric potentials between regions of fixed voltages 
in two dimensions. These programs were modified by B. Eckstein and may now be 
applied to the ion sources and ion collector of a mass spectrometer. Once the 
potentials have been calculated, another computer program produces a histogram 
of ion density versus ion location at a plate within the focusing device. 

From this histogram, the percentage of ion beam transmitted may be found. By 
changing the plate voltages and obtaining a histogram for these voltages, the 
best set of plate voltage values can be determined for maximum ion transmission 
and focusing. 
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II. Ion Focusing Devices 

Four different ion focusing devices were analyzed in this study, three of 
which were mass spectrometer ion sources and one of which was a mass spectro- 
meter ion collector assembly. The first ion source was developed by A. J. 
Loveless and R. D. Russell and will henceforth be referred to as the Loveless 
ion source. Figure 1 is a graph of a two-dimensional cross-section of the 
Loveless ion source. The second ion source is shown in Figure 2 and was 
developed by the California Institute of Technology. Accordingly, it is 
called the Cal Tech source. The last ion source is the National Bureau of 
Standards, or NBS, ion source depicted in Figure 3. The histograms were most 
often produced for an ion of strontium, specifically ®®Sr. 

The Loveless ion source was developed for maximum transmission of a Pb'*’ 
ion beam (1). The 1( ns system is only 2 cm in length and the total accel- 
erating voltage is 5000 V. The Pb"*^ ions have an energy of approximately 0.15 
eV, with a corresponding velocity of magnitude of 4.58 x 10^ cm/s; and the 
center filament is operated at 1750 K. The focus electrodes are at 4995 V and 
the collimating slits at zero volts. The source was studied to find the best 
operating voltages for emission of ^^Sr. 

The Cal Tech ion source is a device that was partially studied last 
summer as part of an LPI undergraduate internship, when the electric potential 
was computed for a maximum voltage of 10. 0 V. The first drawing out plate was 
at 8.5 V, the second at 3.0 V. The collimator was at 1.5 V and the exit slit 
at zero volts. This summer the intermediate voltages for = 10.0 V were 

changed and plate voltages were also scaled to 'Jmax ” 7500 V. Again, ®®Sr was 
used in all histograms. 
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The NBS ion source is the device used in a mass spectrometer run by 
Dr. L. E. Nyquist in Building 31 at the NASA Johnson Space Center. One of the 
elements most commonly analyzed for isotopic ratios is ^^Sr (which is why ^^Sr 
is used in all histogram studies). The maximum voltage is 7500 V with an 
operating temperature of 1300 to 2000^C. The plate voltages were measured 
last summer while the mas.s spectrometer was in operation with a seven digit 


voltmeter and the values obtained were: 

a. Drawing out plate 7484 V 

b. First discriminator 7285 V 

c. Focus plates 6673 V 

d. Second dincriminator 449 V 

e . Exit slit 0 V 


This source was analyzed last summer for a version scaled to a maximum voltage 
of 1000 V, and plots of the equipotential contours and exit slit histograms 
were obtained. This summer the electric potentials were calculated for the 
working voltages and histograms of the individual plates obtained. The 
voltages were varied and the appropriate histogram examined to determine 
whether or not focusing was improved. 

The ion collector assembly consists of a series of plates terminating in 
a photomultiplier tube and is the collecting device for the mass spectrometer 
used by Dr. Nyquist (sec Figure 4). The collector was analyzed to determine 
whether or not a significant percentage of particles were lost before they 
reached the photomultiplier tube and, if so, how to pervent this loss. 
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III. Experimental Procedures 

In order to analyze any focusing device, all dimensions of the lent 
system must first be found, either from blueprints of the device or by actual 
measurement. Next, one must choose which two-dimensional cross-section to 
study. For this project, the cross-section with the smallest apertures was 
chosen because of the more stringent focusing requirements. Then the two- 
dimensional cross-section is placed on graph paper as in Figures 1 through 4. 

The computer programs used in calculating the electric potential and 
histograms are written so that every focusing device must have the dimensions 
of IOC X 100. Thus, once the device has been placed on graph paper, the X and 
Y coordinates must be changed from centimeters (or inches) to integer steps 
from 1 to 100 (see Figures 1-4). From this 100 x 100 array, a boundary 
conditions file may be created that contains the locations of each plate 
voltage and the value of that voltage. File 1 is an example of one Loveless 
ion source. The boundary conditions file is called NAME.BND, where NAME 
refers to the lens system being studied. The suffix BND remains the same no 
matter what device is being studied in order to indicate that this is a 
boundary conditions file. 

The initial velocities of the ions must also be determined. Since 

the ions have an energy of approximately 0.10 eV, the magnitude of the velocity 
is 4.69 X 10^ cm/s. Because the ions may be emitted at any angle with respect 
to the filament, the computer calculates ion trajectories from 0 to 180 degrees 
in 15 degree steps relative to the filament. These X and Y velocities are 
contained in a file called MASS2.VEL (see File 2). Thus, 13 trajectories may 
be computed at any location on the filament. For the ion collector, the 
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angular dispersion of the ^^Sr beam is one degree, so trajectories are 
computed for 89.0 to 9^.0 degrees in 0.2 degree steps at every location from 
which ions originated. The iona have been accelerated through 7300 V and now 
have a velocity of A. 06 x 10^ cm/s. 

Once the boundary conditions file NAME.BND has been created, the inter- 
mediate electric potential of the device may be calculated by means of two 
programs called BNDSET and ELECT. These programs are located in Room 180 of 
Building 31 at the Johnson Space Center on the hard disk numbered RK-29 and 
are designed to run on the building's PDF 11/40 computer located in the same 
room. 

The program BNDSET creates two new files, NAME. RAN, an array file, and 
NAME. BIT, a bitmap file, from the boundary conditions file. (The user may 
call the array and bitmap files anything she or he desires, but retaining the 
boundary conditions file prefix and changing only the suffix reduces confusion.) 
The suffixes RAN and BIT always refer to an array file and a bitmap file, 
respectively. The bitmap contains the location and value of the boundary 
conditions. The array file ultimately will contain the value of the potential 
at each point in the 100 x 100 grid of the device. 

The program BNDSET will also ask if the array file is new or old. If one 
has never calculated the electric potentials for a device of this geometry, 
the response is "N" for new. The computer then initializes each point in the 
array file to zero. If one has already computed the electric potentials for a 
device of this geometry and now desires the electric potentials when a few of 
the plate voltages are changaed, one answers "0" for old. Then the computer 
will correct the old data for the new boundary voltages without having to 


lecalculate Che entire array, and run time is shortened. VIhen Che program has 
finished, the computer prints a dot. DNDSET takes upwards of a minute to run. 

For an actual run of BNDSET, see Program 1. 

l*he actual potential values are produced from Che program called ELECT, 
which may only be started after BNDSET has been run for Che current boundary 
conditions. ELECT asks for Che array file and bitmap file names and asks 
whether Che field is planar or radial (1 or R). Since all focusing devices 
have rectangular apertures, Che fields are planar. ELECT will also ask if one 
wants NAME. RAN printed. Tne usual response is "N" for no, as one's data consists 
of Che histogram. Finally, ELECT will ask if Che current limits of I ■ 2-99 
and J ■ 2-99 are acceptable, which is answered by "Y" (i.e., yes). ELECT sets 
the perimeter potential at zero and calculates the potential for each point 
within the boundary. 

The program ELECT uses a relaxation method Co compute potential values. 

For each pass, Che line printer writes (1) the number of tiroes a pass has been 
made, (2) Che X and Y coordinates of the potential being examined, (3) Che new 
potential value at this location, and (4) the absolute value of Che difference 
between the old potential value and the new. The program only stops calculating 
if one types Cntrl C, or if the difference between Che old and new potential 
values is less Chan 0.001 volts. If one has typed Cntrl C and now wants Co 
continue livaluaCion of the electric potentials so that they may be more 
accurately determined, one need only restart ELECT once morc; and the program 
will Cake up where it left off. One need not run BNDSET again while Che boundary 
conditions remain Che same. Program 2 shows an actual ru” of ELECT. 
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If one is uncertain about the boundary conditions one has chosen and 
wants to be sure that the plates are at their proper locations, one may run a 
program called BMPNT and examine its line printer output. This output consists 
of a series of dots and X's, where Che dots are points in Che bitmap Chat have 
not been set at a voltage and the X's are the locations of the fixed voltages. 
BMPNT only requires the maximum value of IDIM and JDIM, i.e., 100, and the name 
of the bitmap file. Tliis program is also on RK-29. A run of BMPNT is found 
in Program 3. 

Once the electric potentials have been calculated for a lens system, one 
starts a program C0NPL9 that (1) plots equipotential contours, (2) plots ion 
trajectories, and (3) computes a histogram of the Y-location of each particle 
when Che X-coordinaCe reaches a certain value. This program may be found on 
disk RK-26 in room 180 of Building 31 at JSC. (Disk RK-26 also contains all 
boundary conditions, array, and bitmap files for each focusing device.) 

Program C0NPL9 will only run if a plotter is hooked up and is on. It 
first asks for the array file names (see Program A). Then it asks "Contours 
or Trajectories." If one types "C," the computer ask the name of a batch 
input file. A batch input file is simply a list of Che values of Che equipo- 
tential contours to be plotted and is created using Che editor. Aji example of 
a batch input file is found in File 3. All batch input files are located on 
disk RK-26. Once Che contours have been plotted, Che computer again asks 
"Contours or Trajectories." One may then type "C" and reference another batch 
input file, which will result in more contours being drawn on the plot; or one 
may type "T" so that ion trajectories may be calculated. 
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When CrajecCories are to be computed, the computer first asks for the 
bitmap corresponding to the array file already given. It then asks for the 
mass of the ions and a quantity called "cell size." This cell size is the 
length in centimeters that is covered by one unit of the 100 x 100 array in 
the X and Y directions and may be determined from the lens system graph (see 
Figures 1 • 4). Due to a hidden error in the logic of C0NPL9, the mass must 
be entered as a negative number to produce accurate results. 

The computer will next want to know whether or not trajectories are to be 
computed one at a time (Y or N). If the answer is Y, the computer will require 
the initial coordinates XO and YO and the initial velocities VXO and VYO. If 
the answer 's N, the XO and YO values asked for are now starting values. At 
XO and YO, trajectories are computed for ions of any X and Y velocity by 
referencing a file such as MASS2.VEL when the computer asks for the initial 
conditions file. Wlien all ion trajectories in MASS2.VEL have been computed at 
XO and YO, YO is increased by adding to it a quantity called AMT (which C0NPL9 
needs to be told), and trajectories are again computed from MASS2.VEL. This 
process is repeated until the Y-coordinate becomes greater than YUPPER (a 
quantity one must also tell the computer). Then the computer stops calcu- 
lating trajectories. In this way, the computer will literally compute thousands 
of trajectories after being given only five parameters, i.e., XO, YO, YUPPER, 
AMT, and u velocity file. 

C0NPL9 also gives a choice of data output as data may go to the plotter 
and/or to the line printer. If data only goes to the plotter, the only output 
is a plot of each ion trajectory. If data goes to the printer, the computer 
needs additional parameters. First, it requires the change of X for which the 


9 


line printer outputs data. If one answers, for example, "63,83," the line 
printer will print the time, X- and Y-coordinates , and X- and Y-velocities for 
each step of the trajectory when the X*coordinatc is between 63 and 83. If 
one answers "-1," none of the particle's positions are listed. 

If data goes to the line printer and trajectories are computed from a 
velocity file, a histogram is produced once all the ions' paths have been 
calculated. The input "Range of Y for histogram" is the region in the Y 
dimension for which the histogram is computed. One usually chooses the Y* 
range to be symmetric about the center of the lens system. The next input is 
the value of X for which the histogram is computed. This X value is usually 
one array unit above the plate being studied so that one may determine the 
width of the ion beam along with percent transmission. The histogram itself 
is 232 units wide. The first and 232nd locations tell how many particles 
landed beyond the Y-range. The line printer prints the number of particles in 
each unit of the histogram and then gives a graph of the histogram. For an 
actual histogram, see Table 1. 

Once a histogram has been obtained at a plate within the lens system, one 
can calculate the percentage of transmission of the ion be.]m. One may then 
change the voltage on this plate, recalculate the electric potential distri- 
bution, and, by obtaining another histogram, compare transmissions to see 
which is the better voltage. This process may be repeaf=:d until the optimum 
operating voltages have been determined. 
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IV. Results 

A number of factors obstructed complete analyses of the four ion focusing 
devices currently under consideration. First and foremo.<tt was a lack of 
computer time. The set of programs used in this project only run on the POP 
11/40 located in Building 31, and this computer is greatly in demand, especially 
during the summer when special projects requiring computer time on this machine 
are being pursued by LPI interns and NASA personnel. Also, the POP 11/40 is 
designed for running equipment and not for batch work and thus is not as fast 
as a computer specifically designed for these types of programs. 

A third problem lay in the trajectory plotting program itself. Last 
summer the version used wa.? called C0NPL7. This program produced all the data 
that C0NPL9 gives, but did not require the cell size. As a result, the initial 
velocities of the ions were arbitrarily scaled to the array without regard for 
the source's actual dimensions; and all histograms produced by CONPL7 last 
summer are wrong. (The CONPL? programs are derived from a program CONP20 
written by Dan Anderson a few years previously and were an aim of this project, 
i.e., to produce a program capable of evaluating the focusing properties of a 
lens system.) Initially, C0NPL7 was changed to a version called C0NPL8, where 
CONPL8 included the array cell size. In CONPL8, however, it was discovered 
that the step size in the trajectory calculation was too large to produce an 
accurate histogram. Accordingly, the step size was reduced considerably (by a 
factor of 3) so that trajectories of 0.2% accuracy are obtained; and the new 
program was named CONPL9 for obvious reasons. Since C0NPL9 takes much longer 
to calculate a trajectory, the velocity file was decreased from one containing 
61 trajectories to one containing only 13 for the mass spectrometer sources. 
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Another delay in getting data reaulta from the inconsistencies of the 
PDP 11/40. The progiam C0NPL9 requires a line printer for the histogram, and 
the computer in Building 31 is the only 11/40 available with a line printer. 

Thus, when the computer went down, the project had to stop also, because 
almost all work is done by computer. (The histograms take very little time to 
analyze.) The only person capable of fixing the PDP 11/40 on short notice 
when it dies is Dan Anderson, who was on vacation for three weeks this suimner. 

While he was gone, the PDP 11/40 passed away twice, and over a week of work 
time was lost. 

However, a few results were obtained for the lens systems studied. For 
each source, a histogram at the plate whose voltage was being changed was used 
Co determine the best plate voltage. However, as indicated by Figures 5 and 
6, Che criterion should have been Che transmission through Che exit slit, 
regardless of the plate being changed. 

For each focusing device, a plot of the equipoCential contours corres- 
ponding Co Che initial voltages was obtained. Then a few trajectories were 
plotted. Finally, the histograms were obtained for Che different voltages (see 
Figures 5-8). 

For the Loveless ion source, one can see from Figures 5 and 6 that, 
although 5000 V on the focusing electrodes gives 100% transmission through the 
electrodes, no ions pass through the exit slit. The data charts for each 
histogram are shown in Tables 2 and 3. 
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The C«1 Tech source has hardly been analyzed at a'l but the data in 
Table U seems to indicate that the transmission through the first drawing out 
plate will be improved if the voltage is increased further. 

The NBS ion source apparently requires the least improvement, although an 
exit slit histogram was not obtained for the new drawing out plate voltage cf 
7498 V (see Table 5). However, without changing the voltages on the first 
discriminating plate and focus plates, the entire beam transmitted by the 
drawing out plate is passed through the focus plates, i.e., 73. OZ of the total 
ion beam is transmitted through the focus plates. 

The ion collector assembly has been analyzed the least because of its 
efficiency, as it delivers 84.8% of the total ion beam to the conversion dynode 
of the photomultiplier tube. Since the collector assembly has e lower priority 
than the ocher focusing devices, it was not analyzed further. 

At present, exit slit histograms are being computed for each of the three 
ion sources being studied. The data, however, has not yet been analyzed. 
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V. Further Possibilities 

As one can tell by the (so far) sketchy results of this project, there is 
much more work to be done. First and most importantly, all four focusing 
devices must be completely analyzed by getting a histogram for the exit slit 
after the new electric potential distribution has been calculated. 

The experimental procedure may be improved in a number of ways. One 
possible improvement is to put the programs onto different computer systems at 
JSC and LPI so that different lens systems may be analyzed simultaneously. 
Also, the user is no longer dependent upon the vagaries of the POP 11/40. The 
project would also benefit if the programs were placed on faster computers, 
such as the Univac system in Building 12. 

One may also want to determine whether or not it is possible to combine 
the three programs BNDSET, ELECT, and CONPL9 into a single program that will 
automatically change plate voltages until the best focusing conditions are 
obtained. A computer faster than the 11/40 would be a necessity, however, as 
the combined programs would take weeks to run otherwise. 

Another variation of an ion source that has yet to be analyzed is one 
that uses a V-shaped filament. The first draft of the commands necessary to 
change C0NPL9 so that a V-shaped filament may be analyzed was written this 
summer but was never tested. 

The ion collector assembly should also be analyzed for the electron 
trajectories from the conversion dynode of the photomultiplier tube. 

Any of these ion sources may also be analyzed in three dimensions from 
programs written by Dan Anderson. These programs were not used, this summer 
because heir run time is longer. 

Finally, one may experiment with entirely new lens systems that have 
different geometries from any focusing devices currently in use. 



VI. Conclusions 

The data for each of the four focusing devices is obviously incomplete 
and can be improved if the computer time is available. One interesting result 
is that exit slit transmission is not necessarily improved by increasing 
individual plate transmission. Thus, one should always look at the exit slit 
histogram instead of the histogram for the plate whose voltage is being c’janged, 
otherwise, the data may be misleading. 
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